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METHOD AND SYSTEM FOR 
SEGMENTING MAGNETIC 
RESONANCE IMAGES 

Background of Invention 

[0001] The invention relates to magnetic resonance imaging (MRI) and image processing 
methods. More particularly, the invention relates to segmenting brain images into 
various structures found within the brain. 

[0002] It is well known to employ MRI for visualization of internal anatomical structures 
such as brain imaging. Brain imaging, and particularly for use in brain volume 
measurement is of particular interest for studying degenerative brain diseases such as 
rf Alzheimer's disease (AD). Examination of the Alzheimer brain pathology shows 

1'U extensive p -amyloid plaque, neuron tangles and brain atrophy. Typically, magnetic 

g resonance imaging brain volume measurements are used to monitor the disease 

^ progression. Normal aging brain atrophy is only about a 3.5% decrease per decade, 

but the rate of atrophy increases in subjects exhibiting dementia. Thus, brain volume 
measurements provide a measurable correlation available to predict Alzheimer T s 
disease. The diagnosis of other brain diseases such as tumors and edema also rely on 
brain volume measurement. 

[0003] 

Segmentation techniques typically require some user interaction such as placing a 
seed point within an area of interest in order to initiate segmentation or registration 
processing. Further, in order to segment multiple structure or multiple tissue types 
within the brain generally requires repeating the seed placement and subsequent 
processing for the selected structure or tissue type. Thus, measuring brain volume 
from magnetic resonance images either by registration methods or by segmentation 
methods is typically tedious because there is considerable manual editing and 
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processing needed. 

[0004] Measurement of brain structures could lead to important diagnostic information 
and could also indicate the success or failure of a certain pharmaceutical drug. A 
necessary step to quantify brain structures is the availability of a segmentation 
technique of some sort. Other techniques, such as neural network techniques require 
a training stage for the algorithm. These techniques are time consuming, which make 
th em less clinically desirable. Therefore, there is a need for a fast and accurate 
unsupervised segmentation technique. 

[0005] What is needed is a method and system for segmenting medical images in an 

automatic manner with minimal user intervention. What is further needed is a method 
□ for segmenting medical images for a plurality of structures within an anatomical body 

^ of interest. 

J Summary of Invention 

m [0006] In a first aspect, a method is provided for segmenting magnetic resonance (MR) 
~ images of an anatomical body of interest, for example the brain. The method 

Q comprises classifying a plurality of selected structures within the body of interest 

I y based on a plurality of image processing computations wherein the computations are 

^ adapted to relate respective T2 relaxation times corresponding to each of the 

fy structures. Thereafter, the method comprises segmenting the MR images for each of 

the structures substantially concurrently based on the plurality of image 

computations. 

[0007] In a second aspect, a system is provided for automatically segmenting magnetic 
resonance (MR) images of an anatomical body of interest. The system comprises a 
processor coupled to an MR image acquisition device wherein the processor is 
adapted to perform concurrent segmentation computations for a plurality of selected 
structures within the anatomical body of interest. The system further comprises an 
interface unit coupled to the processor adapted to present information relating to the 
segmented computations corresponding to the plurality of selected structures. 

Brief Description of Drawings 

[0008] The features and advantages of the present invention will become apparent from 
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the following detailed description of the invention when read with the accompanying 
drawings in which: 

[0009] Figure 1 is a block diagram illustration of a medical imaging system for which 
embodiments of the present invention are applicable; 

[001 0] Figure 2 is a schematic illustration of an exemplary embodiment of a method for 
segmenting MR images; 

[0011] Figures 3 and 4 are dual echo images of types useful in the method of segmenting 
MR images shown in Figure 2; 

[001 2] Figures 5 and 6 are graphical illustrations of a scatter plot and radial histogram, 
respectively, that are generated from the images of Figures 3 and 4 and useful in the 
method of segmenting MR images shown in Figure 2; and, 

[001 3] Figures 7-1 2 are graphical illustrations of segmented brain structures and their 

corresponding T2 windows derived by the method of segmenting MR images shown in 
Figure 2. 

Detailed Description 

[001 4] Referring to Figure 1 , a general block diagram of a system 1 00 for disease 
detection is shown. System 100 includes an imaging device 110, which for the 
embodiments of the present invention is a magnetic resonance imaging (MRI) system 
used to generate a plurality of MR images. 

f001 S] During a MR imaging session, the patient is placed inside a strong magnetic field 
generated by a large magnet. Magnetized protons within the patient, such as 
hydrogen atoms, align with the magnetic field produced by the magnet. A particular 
slice of the patient is exposed to radio waves that create an oscillating magnetic field 
perpendicular to the main magnetic field. The slices can be taken in any plane chosen 
by the physician or technician (hereinafter the "operator") performing the imaging 
session. The protons in the patient's body first absorb the radio waves and then emit 
the waves by moving out of alignment with the field. As the protons return to their 
original state (before excitation), diagnostic images based upon the waves emitted by 
the patienf's body are created. The MR image slices are reconstructed to provide an 
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overall picture of the body area of interest. Parts of the body that produce a high 
signal are displayed as white in an MR image, while those with the lowest signals are 
displayed as black. Other body parts that have varying signal intensities between high 
and low are displayed as some shade of gray. 

[001 6] Once initial MR images have been obtained, the images are generally segmented. 
The segmentation process classifies the pixels or voxels of an image into a certain 
number of classes that are homogeneous with respect to some characteristic (i.e. 
intensity, texture, etc.). For example, in a segmented image of the brain, the material 
of the brain can be categorized into three classes: gray matter, white matter, and 
cerebrospinal fluid. Individual colors can be used to mark regions of each class after 
the segmentation has been completed. Once the segmented image is developed, 
surgeons can use the segmented images to plan surgical techniques. 

[001 7] F|C 1 j|| ustrates a me dical imaging system 100 to which embodiments of the 

invention are applicable. The system includes an imaging device 1 10, a processor 120 
and an interface unit 130. Imaging device 110 is adapted to generate a plurality of 
image data sets 140 and is, for example, a magnetic resonance (MR) scanner. In the 
context of MR, acquisition of image data is generally referred to as "scans". Processor 
120 is configured to perform computations in accordance with embodiments of the 
present invention that will be described in greater detail with reference to Figures 2- 
1 2. Processor 120 is also configured to perform computation and control functions for 
well-known image processing techniques such as reconstruction, image data memory 
storage, segmentation and the like. Processor 120 may comprise a central processing 
unit (CPU) such as a single integrated circuit, such as a microprocessor, or may 
comprise any suitable number of integrated circuit devices and/or circuit boards 
working in cooperation to accomplish the functions of a central processing unit. 
Processor 1 20 desirably includes memory. Memory within processor 1 20 may 
comprise any type of memory known to those skilled in the art. This includes Dynamic 
Random Access Memory (DRAM), Static RAM (SRAM), flash memory, cache memory, 
etc. While not explicitly shown in FIG. 1 , the memory may be a single type of memory 
component or may be composed of many different types of memory components. 
Processor 1 20 is also capable of executing the programs contained in memory and 
acting in response to those programs or other activities that may occur in the course 
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of image acquisition and image viewing. As used herein, "adapted to", "configured" 
and the like refer to mechanical or structural connections between elements to allow 
the elements to cooperate to provide a described effect; these terms also refer to 
operation capabilities of electrical elements such as analog or digital computers or 
application specific devices (such as an application specific integrated circuit (ASIC)) 
that are programmed to perform a sequel to provide an output in response to given 
input signals. 

[001 8] Interface unit 1 30 is coupled to processor 1 20 and is adapted to allow human 

users to communicate with system 1 00. Processor 1 20 is further adapted to perform 
computations that are transmitted to interface unit 1 30 in a coherent manner such 
that a human user is capable of interpreting the transmitted information. Transmitted 
information may include images in 2D or 3D, color and gray scale images, and text 
messages regarding diagnosis and detection information. Interface unit 130 may be a 
personal computer, an image workstation, a hand held image display unit or any 
convention image display platform generally grouped as part of a CT or MRI system. 
Referring further to Figure 1 , processor 1 20 is also adapted to perform the 
embodiments of segmentation methods that will be described in greater detail with 
reference to Figures 2-1 2. 

[0019] All data gathered from multiple scans of the patient is to be considered one data 
set. Each data set can be broken up into smaller units, either pixels or voxels. When 
the data set is two-dimensional, the image is made up of units called pixels. A pixel is 
a point in two-dimensional space that can be referenced using two-dimensional 
coordinates, usually x and y. Each pixel in an image is surrounded by eight other 
pixels, the nine pixels forming a three-by-three square. These eight other pixels, 
which surround the center pixel, are considered the eight-connected neighbors of the 
center pixel. When the data set is three-dimensional, the image is displayed in units 
called voxels. A voxel is a point in three-dimensional space that can be referenced 
using three-dimensional coordinates, usually x, y and z. Each voxel is surrounded by 
twenty-six other voxels. These twenty-six voxels can be considered the twenty-six 
connected neighbors of the original voxel. 

[0020] 

Referring to Figure 2, an embodiment for segmenting MR images is provided that 
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segments and quantifies brain structures from T2 dual echo MR images. As used 
herein, n T2", "T2 parameter" and the like refer to the time constant, or alternatively 
spin-spin relaxation time, T2 that is well known in the art of MR imaging. T2 is the 
time measurement for a given nuclei to return to be uniformly distributed around the 
static magnetic field (referred to as "B") once the RF pulse sequence is completed in 
the MR scan. There is a T2 value associated with a given tissue type or brain structure, 
thus the T2 value is useful in distinguishing selected tissue types in a MR image. 
Further, the given 12 value may be visualized differently between dual echo images. 
For example, the CSF typically has higher values in the second echo and the face has 
higher values in the first echo. In an exemplary embodiment, four different tissue 
clusters are classified in a scatter plot (air, cerebrospinal fluid (CSF), brain, and face). 
In further embodiments, additional image-processing techniques were implemented 
to reduce the spread of these clusters and subsequently generate tissue based 
information datasets (hereinafter referred to as "T2 windows") that encompass all the 
information needed to segment and subsequently quantify the corresponding tissues 
in an automatic fashion. 

[0021 ] The input to the method shown in Figure 2 are images acquired at step 21 0 by 

MRI scanning, for example on a 1 .5 Tesla CE Signa MR scanner. The MR data consists 
of dual echo coronal images of the head. The dual echo was acquired by known 
methods using T2 spin echo pulse sequence. In an exemplary embodiment, the first 
echo, proton density weighted (PDW), was acquired at 30 msec and the second echo, 
T2 weighted (T2W), was acquired at 80 msec. The voxel size was .9375 x .9375 x 4 
mm3. Figure 3 illustrates the PDW image and Figure 4 illustrates the T2W image. It is 
to be appreciated by those skilled in the art that other dual echo pulse sequences may 
also applicable to methods described herein. 

[0022] Before describing the segmentation method of the present invention as shown in 
Figure 2, a brief description of the image processing steps will be described in greater 
detail below. Each of the image processing filters is applied to dual echo images 
shown in Figures 3 and 4 as described above. Referring to Figures 5-6, the main 
image processing steps are a scatter plot (SP), a T2 radial histogram (RH), a T2 median 
filter, and a dual echo MIP filter. Each of these steps is defined herein as an image 
processing computation based on T2 values and other information relating to T2. 
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[0023] Referring first to Figure 5, the SP is a two-dimensional (2D) histogram constructed 
using the dual echo MR images. The abscissa consists of the range of voxel values of 
the first echo and the ordinate consists of the range of voxel values from of the 
second echo. The intensity value of a pixel in the scatter plot represents the frequency 
of occurrence of voxel values in the PWD image that are related to voxel values in T2W 
image at corresponding locations. 

[0024] Referring to Figure 6, the T2-RH is a histogram plotting the linear profile angle 
versus the area under the profile. The abscissa in the T2-RH represents the angle in 
the scatter plot. For instance, for zero-degree angle the profile follows on the 
abscissa of the SP. For 90-degree angle the profile follows the ordinate of the SP. The 
ordinate in T2-RH represents the sum of the intensity values (area under the profile) 
Q for a linear profile. The profile starts from the lower left in the scatter plot with an 

^ angle to the edge of the scatter plot. Note that the first 20 values in the profile were 

~ skipped because they represent part of the air cluster. A linear profile in the SP 

gi represents all the pixels in the corresponding brain image that have the same T2 

;jf values. T2 values extend along the profile mainly because of the difference in 

composition of proton density between tissues. T2 values for a tissue also extend in 
(2 an angular fashion, thus forming clusters, mainly because of MR signal 

j Jf inhomogeneity, volume averaging, and coil shading artifacts. 

f[j [0025] Referring further to Figure 2, the T2 median filter step 220 is first applied to the 
dual echo images and is adapted to reduce the spread of the T2 relationship values 
between the dual echo images. The T2 Median Filter applies the median of the 
calculated T2 values from the dual echoes instead of using the median of the pixel 
values from either echoes. The inputs to this filter are dual echo images (PDW and 
T2W images) and the outputs are dual echo filtered images that retain the MR 72 
relationship. Note that the output of this step is a dual echo filtered data that 
maintains the tissue MR T2 parameter. 

[0026] 

The dual echo maximum intensity projection (MIP) steps 230 and 240 adapted to 
enhance and amplify one peak distribution while suppressing other peaks in T2-RH. 
Further, the dual echo MIP is adapted to maintain the spatial relationship between the 
dual echo images. In an embodiment for the dual echo MIP, processing selects a 
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preferred echo and implements a MIP on it and then provides the spatial location of 
the implemented MIP. The spatial location is then used to extract values from the 
other echoes. Thus, each voxel retains the T2 characteristics between the MR echoes. 
The dual echo MIP can be based on the PDW image (230) or the T2W image (240). The 
dual echo PWD-MIP is adapted to enhance the face peak distribution while 
suppressing the cerebrospinal fluid (CSF) and the brain distributions in the T2-RH. 
Note that this filter outputs two echoes. In an exemplary embodiment, the first echo is 
the conventional MIP of the PWD images. The second MIP echo is generated using the 
MIP pixel locations from the PWD image to get the corresponding pixel values from 
the T2W image. The dual echo T2W-MIP outputs the second echo with the 
conventional MIP, while the first echo is generated using the MIP pixel location from 
the T2W image to get the corresponding pixel values from the PDW image. The MR T2 
parameter is maintained in the dual echo MIP output images. 

[0027] An embodiment for a method for segmenting magnetic resonance (MR) images of 
an anatomical body of interest, for example the brain, comprises classifying a plurality 
of selected structures within the body of interest based on a plurality of image 
processing computations as described above wherein the computations are adapted 
relate respective T2 relaxation times corresponding to each of the structures. 
Thereafter, the method comprises segmenting the MR images for each of the 
structures substantially concurrently based on the plurality of image computations. In 
an embodiment described herein, the anatomical body of interest is the brain and the 
plurality of structures include at least one of air, face tissue, brain tissue, 
cerebrospinal fluid (CSF), edema and tumor tissue. Desirably, the segmenting step is 
used in at least one of diagnosis, volume measurement and clinical research. As 
described above, the plurality of image processing computations comprises a scatter 
plot of voxel values of the MR images, at least one radial histogram, and a plurality of 
image processing filters. The image processing filters are adapted to enhance T2 
values of the MR images and further relate dual echo image information. 

[0028] T2 median filter step 220 is adapted to reduce the spread of the T2 relationship 
values between the dual echo images, as described above. The output of the T2 
median filter is a dual echo filtered data set that is used as an input to the dual echo 
PWD-MIP step 230 as well as the T2W-MIP step 240. Next, three scatter plots (SP0, 
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SP1 , and SP2 shown as steps 250, 260 and 270, respectively) are generated and three 
T2 radial histograms (RHO, RH1 , and RH2 shown as steps 280, 290 and 300, 
respectively in Figure 2). SPO and RHO are generated directly from the filtered dual 
echo data. SP1 and RH1 are generated from the output of the dual echo PWD-MIP. SP2 
and RH2 are generated from the output of the dual echo T2W-MIP. 

[0029] The scatter plot and radial histogram information is used to detect peaks for the 
various brain structures. Using the gradient descent technique, the peak of the brain 
distribution (Dl) is detected from the RHO, the peak of the face distribution (DO) from 
RM, and the peak of the CSF distribution (D2) from RH2. From the angular profiles 
DO, Dl, and D2, the peaks of these profiles are also determined using also a gradient 
descent technique. These peaks correspond to the centers of the tissue clusters (face, 
brain, and CSF) in the scatter plot, as shown in Figure 5. The center for the air 
distribution is located in the lower left of the scatter plot. The tissue clusters 
described above are used to generate an information dataset for use in segmenting 
the various tissue types. 

[0030] The output resulting from processing steps 210-300 shown in Figure 2 results in 
an information dataset (hereinafter "T2 window") for each selected tissue type or 
structure. The information dataset, T2 window, is based on the image processing 
filters described above. The image processing filters are each adapted to extract, or 
alternatively enhance, various aspects of the T2 values corresponding to the selected 
structures of the brain. Referring to Figures 7-1 2, there is shown a graphical 
illustration of segmented brain structures and their corresponding T2 window. Figures 
7 and 8 illustrate the segmented CSF and CSFT2 window, respectively. Figures 9 and 
10 illustrate the segmented brain and brain T2 window, respectively. Figures 1 1 and 
12 illustrate the segmented face and face T2 window, respectively. 

[0031] 

Each of the T2 windows appears like an annulus region of interest that is 
determined by four parameters. These are the lower and upper angular limits, and the 
inner and outer circles. In an exemplary embodiment, the lower angular limit for the 
CSF T2 window was determined by averaging Dl and DO, where the upper limit was 
set to 60 degrees to compensate for noise and artifacts. Note that the maximum angle 
for any point with information in the scatter plot is expected to be at 45 degrees since 
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the intensity values in the T2W image are always lower than the corresponding 
intensity values in the PDW image (MR T2 relaxation). The inner circle was set with a 
radius equal to half the distance between the center of the CSF cluster and the center 
of the air cluster. The outer circle was set to be three times larger than the inner 
circle. The face and brain T2 windows were determined with the same technique. 

[0032] The CSF was segmented and quantified using the T2 windows from all of the 

corresponding voxels in the original MR images (Figures 7-1 2). However, the brain T2 
window was used to segment the corresponding voxels from the original MR images. 
In further embodiments and with further reference to Figure 2, desirably an erosion 
filter 330, then region growing 340 with a seed taken from the center brain cluster, 
and finally dilation 350 are applied to segment the brain images. These additional 
, sa% steps remove misclassified voxels (typically the lower face area) that have similar brain 

5 *P T2 values. Subsequently, the volume was calculated for the brain and face. 

: s ;J [0033] In embodiments of the present invention, the T2 window datasets provide the 
CM needed information to segment the CSF, brain, and face in a totally automatic, and 

j ; I alternatively unsupervised, manner. Further, the T2 window computations for each 

structure within the brain are computed simultaneously, or alternatively concurrently, 
\*A therefore enabling presentation of each of the segmented structures after processing. 

; if It is to be appreciated that the segmented structures are processed concurrently, thus 

□ reducing image processing time. Compared with conventional segmentation 

techniques requiring seed placement in a given structure and then waiting for 
subsequent segmentation processing, embodiments of the method shown in Figure 2 
provides the ability to segment a plurality of structures or tissues types automatically 
and substantially simultaneously. 

[0034] Embodiments of the present invention have been described in the context of 

diagnosis and brain volume measurement. Further embodiments would be applicable 
to clinical research applications where there is a need to process multiple images and 
segment the multiple images for various structures. By employing the image 
processing computations described above, the tedious and manual intervention 
generally employed clinical research of multiple images is substantially reduced. 

[0035] Further embodiments of the image processing computation would be adapting the 
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computations to segment other structures such as edema, tumors in brain imaging. 
The embodiments described herein segments only the four mentioned structures (air, 
face, brain, and CSF). For instance, it does not differentiate between white matter from 
gray matter. The dual echo pixel values for these structures fall within the brain T2 
window because the T2 difference between these structures is too small. Furthermore, 
if the brain has other structures such as edema of because of present head injury or 
tumor, the segmentation methods can be adapted to include the a special 12 edema 
window to segment edema. Similarly, the segmentation methods can be adapted to 
include a T2 window to segment for brain tumors. In further embodiments, the 
method for segmenting is adaptable to other anatomical bodies (e.g. heart or other 
organ) of interest that may benefit from image processing computations based on T2 
values as described herein. 

[0036] In an embodiment for the system shown in Figure 1 , a system is provided for 

automatically segmenting magnetic resonance (MR) images of an anatomical body of 
interest. The system comprises processor 1 20 coupled to MR image acquisition device 
1 10 wherein processor 1 20 is adapted to perform concurrent segmentation 
computations for a plurality of selected structures within the anatomical body of 
interest. The system further comprises interface unit 1 30 coupled to the processor 
adapted to present information relating to the segmented computations 
corresponding to the plurality of selected structures. 

[0037] While the preferred embodiments of the present invention have been shown and 
described herein, it will be obvious that such embodiments are provided by way of 
example only. Numerous variations, changes and substitutions will occur to those of 
skill in the art without departing from the invention herein. Accordingly, it is intended 
that the invention be limited only by the spirit and scope of the appended claims. 
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